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Abstract 

A unit-vector field n : P — > S 2 on a convex polyhedron P C M 3 
satisfies tangent boundary conditions if, on each face of P, n takes 
values tangent to that face. Tangent unit-vector fields are necessarily 
discontinuous at the vertices of P. We consider fields which are con- 
tinuous elsewhere. We derive a lower bound Ep(h) for the infimum 
Dirichlet energy E l p^{h) for such tangent unit-vector fields of arbitrary 
homotopy type h. Ep(h) is expressed as a weighted sum of minimal 
connections, one for each sector of a natural partition of S 2 induced by 
P. For P a rectangular prism, we derive an upper bound for E^^h) 
whose ratio to the lower bound may be bounded independently of h. 
The problem is motivated by models of nematic liquid crystals in poly- 
hedral geometries. Our results improve and extend several previous 
results. 
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1 Introduction 



S 2 -valued harmonic maps on three-dimensional domains with holes were 
studied in a well known paper by Brezis, Coron and Lieb jlj. As a sim- 
ple representative example, consider the domain Q = R 3 — {r 1 , . . . , r n } (for 
which the holes are points), and let n : Q — > S 2 denote a unit- vector field on 
Q. For Vn square-integrable, we define the Dirichlet energy of n to be 

E(n) = [ {VnfdV. (1.1) 
Jn 

Continuous unit- vector fields on Q may be classified up to homotopy by their 
degrees, d = (d l , . . . , d n ) G Z n , on spheres about each of the excluded points 
(the restriction of n to such a sphere may be regarded as a map from S 2 into 
itself). In order that n have finite energy, we must have that 



0. (1.21 



Let Cn{d) denote the homotopy class of continuous unit- vector fields with 
degrees d satisfying ()1.2|) . and let HCa(d) denote the elements of C^d) with 
finite Dirichlet energy. Let E^(d) denote the infimum of the energy over 
HCn(d), 

E£ i (d)= mf £(n). (1.3) 

It turns out that E^ { (d) is just 8-7T times the length of a minimal con- 
nection on Q. We recall the definition of a minimal connection. Given two 
m-tuples of points in R 3 , P = (a 1 ,..., a m ) and J\f = (b 1 , . . . , b m ) (whose 
points need not be distinct), a connection is a pairing (a- 7 , b 71 ^) of points in 
V and Af, specified here in terms of a permutation n G S m (S m denotes the 
symmetric group). The length of a connection is the sum of the distances 
between the paired points, and a minimal connection is a connection with 
minimum length. Let 

m 

L{V,Af) = min V \ci j - b" (i) | (1.4) 

denote the length of a minimal connection between V and TV. Let \d\ = 
lEil^'l- Then 

Theorem 1.1. j3] The infimum E 1 ^ (d) of the Dirichlet energy of continuous 
unit-vector fields on the domain Q = M. 3 — {r 1 , . . . , r"} of given degrees d? 
about the excluded points r 7 is given by 

E™ f (d) = 87fL{V{d),AT{d)), (1.5) 
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where V(d) is the \d\-tuple of excluded points of positive degree, with r J in- 
cluded d^ times, andN(d) is the \d\-tuple of excluded points of negative degree, 
with r k included \d k \ times. 



In this paper we consider a natural variant of this problem which emerges 
from a boundary- value problem of some physical and technological interest; 
the domain is taken to be a polyhedron on which n is required to satisfy 
tangent boundary conditions. Let P denote a convex bounded polyhedron in 
M 3 , including the interior of the polyhedron but excluding its vertices. Let 
n : P — > S 2 denote a unit-vector field on P. We say that n satisfies tangent 
boundary conditions, or is tangent, if, on each face of P, n takes values 
tangent to that face. (It is clear that this condition could not be satisfied at 
the vertices of the polyhedron, which belong to three or more faces). 

One motivation for the problem comes from liquid crystals applications, 
in which n describes the mean local orientation of a nematic liquid crystal, 
and the Dirichlet energy (jl.6|) coincides with the elastic or Frank-Oseen en- 
ergy in the so-called one-constant approximation (see, eg, |HJ ESI O 122])- 
Polyhedral cells have been proposed as a mechanism for engendering bista- 
bility - they may support two nematic configurations with distinct optical 
properties, both of which are local minima of the elastic energy [31201 QUI- in 
many cases of interest the orientation at interfaces is well described by tan- 
gent boundary conditions, and low-energy local minimisers appear to have 
different topologies. We also remark that harmonic maps between Rieman- 
nian polyhedra have been studied by Gromov and Shoen [2| and Eells and 
Flugede jB] , in particular in cases where the target manifold has nonpositive 
curvature. 

Here we will restrict our attention to continuous tangent unit-vector fields 
on P. (Let us note that, while nematic orientation is, in general, described 
by a director, or PP 2 -valued field, a continuous director field on a simply 
connected domain such as P can be lifted to a continuous unit- vector field.) 
Continuous tangent unit-vector fields on P can be partitioned into homo- 
topy classes Cp(h) labeled by a complete set of homotopy invariants denoted 
collectively by h. A full account of this classification is given in [21] (see 
also |25J). Below we reprise the results we need here. General discussions of 
topological defects in liquid crystals are given in [THl HU H2j • 

For Vn square integrable on P, let 



denote its Dirichlet energy. Let TlCp(h) denote the elements of Cp(h) with 
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finite Dirichlet energy, and let 

E-f(h)= mf £(n) (1.7) 

n&HC P {h) 

denote the infimum of the energy over 7iCp(h). Our first result ( Theorem II. 21 
below) is a lower bound for Ep f (h). This is expressed is terms of certain 
homotopy invariants called wrapping numbers, which we now define. Let / 
denote the number of faces of P, F c the cth face of P, and F c the outward 
normal on F c , where 1 < c < /. For each face we consider the great circle 
on S 2 containing the unit-vectors s tangent to it, ie F c • s = 0. These / (not 
necessarily distinct) great circles partition S 2 into open spherical polygons, 
which we call sectors. The sectors are characterised by sgn (F c • s). We write 

S° = {sG S 2 | sgn (s ■ F c ) =<x c }, (1.8) 

where a = (cr 1 , . . . , a?) is an /-tuple of signs. It should be noted that most of 
the S a, s are empty; indeed, reckoning based on Euler's formula for polygonal 
partitions of the sphere ( (vertices) - |edges| + | faces | = 2) shows that there 
are at most (and generically) f 2 — f + 2 nonempty sectors. 

Next, let C a denote a smooth surface in P which separates the vertex v a 
from the others. For definiteness, take C a to be oriented so that v a lies on 
the positive side of C a . We call C a a cleaved surface. Given n G Cp(h), let 
n a denote its restriction to C a . The wrapping number w aa is the number 
of times n a covers S a , counted with orientation. For n differentiate, this is 
given by 

±;J n'tfu,), (1.9) 

where u is the area two- form on S 2 , normalised to have integral Air, \ a is 
the characteristic function of S a C S 2 , and A a = \J g2 x a w\ is the area of S a . 
Alternatively, w acr can be expressed as the index of a regular value s G S a of 
n a , ie 

w ™ = ^ sgn detrfn a (r). (1.10) 

r | n a (r)=s 

One can show that w au does not depend on the choice of s G S a nor on the 
choice of cleaving surface C a , that the definition (|1.9J) can be extended to 
continuous n, and that its value depends only on the homotopy class of n 
|21j . In fact, the wrapping numbers constitute a complete set of invariants, as 
is shown in the Appendix. They are not independent, however. For example, 
continuity in the interior of P (absence of singularities) implies that, for all 

(/•"" (I. (1.11) 
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where the sum is taken over vertices v a . Continuity on the faces and edges of 
P implies additional constraints. We will say that h = {w acr } is an admissible 
topology if it can be realised by some continuous configuration n : P — > S 2 . 

Theorem 1.2. Let h = {w aa } be an admissible topology for continuous tan- 
gent unit-vector fields on a polyhedron P. Then 

E™ f (h) > Ep(h) :=^2A°L{V a {h),N a {h)), (1.12) 

(7 

where V a (resp. Af a ) contains the vertices of P for which w acr is positive 
(resp. negative), each such vertex included with multiplicity \w acr \. 

Thus, to each sector a may be associated a constellation of point defects 
at the vertices v a of degrees w aa '. The lower bound Ep(h) is a sum of the 
lengths of minimal connections for these constellations weighted by the areas 
of the sectors. 

Theorem ll.2l is proved using arguments similar to those used to show that 
E^ f (d) > 8irL(V(d),Af(d)) in the proof of Theorem ll.il In Theorem II .![ one 
obtains an equality for E^ { (d), rather than just a lower bound, by construct- 
ing a sequence whose energies approach 8irL(V(d),N(d)). It can be 
shown that a subsequence n( fc ) approaches a constant away from a minimal 
connection while |Vn( fe )| 2 approaches a singular measure supported on the 
minimal connection [IJ. In the present case, tangent boundary conditions 
preclude such a construction; n is required to vary across the faces and, 
therefore, throughout the interior of P. However, for P a rectangular prism, 
we can show that Ep(h) correctly describes the dependence of the infimum 
energy on homotopy type. 

Theorem 1.3. Let P denote a rectangular prism with sides of length L x > 
L y > L z and largest aspect ratio k = L x /L z . Then 

Ep f (h) <CK 3 E p (h) (1.13) 

for some constant C independent of h and L x , L y , L z . 

The upper bound of Theorem II. HI is obtained by estimating the energy 
of explicitly constructed tangent unit-vector fields which satisfy the Euler- 
Lagrange equations near each vertex. 

The general form of the Frank-Oseen energy is given by El 122] 

E FO (n) = J [ifi(divn) 2 + K 2 (n ■ curln) 2 + K 3 (n x curln) 2 

+ J FT 4 div((n- V)n - (divn)n)] dV, (1.14) 
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where the elastic constants Kj are material-dependent. It is easily shown that 
tangent boundary conditions imply that the contribution from the K^-term 
in (jl.l4|) vanishes. The elastic constants K%, Ki and K3 are constrained 
to be nonnegative, and the one-constant approximation (jl.fij) follows from 
taking K\ — K 2 = K% = 1. We remark that Theorems 11.21 and 11.31 imply the 
following bounds for the Frank-Oseen energy: 



where K_ (resp. K + ) is the smallest (resp. largest) of the elastic constants 
Kx, K 2 and K 3 . 

Theorems 11.21 and II. HI improve and extend several earlier results. In [T2| 
we obtained a lower bound for E l p { (h) of the form 



where the maximum is taken over f a 's such that |f° - < |v a - v 6 |. The 
quantity (jl.lfij) is generally less than the lower bound given by Theorem II. 2[ 
in particular because it allows for cancellations between wrapping numbers 
of opposite sign. For example, for a regular tetrahedron with sides of unit 
length, ()1.16j) gives a lower bound of Y2aa A a w aa ', whereas Theorem II . 21 gives 
the lower bound ^ a \w aa \. For a rectangular prism, Theorem 11.31 does 
not hold if fTT6j) is substituted for Ep(h). (£□!§ can be directly compared 
to (|2.15j) below, which gives an equivalent (dual) expression for Ep(h).) 
A restricted example of Theorem 11.21 was given in [IBj for the case h is 
a reflection-symmetric topology. These are the topologies of configurations 
which are invariant under reflections through the midplanes of the prism. 

Results related to Theorem 11.31 were obtained for the special case of 
reflection-symmetric topologies in PBJ EE] • The constructions and estimates 
are simpler in this case, and one can show that the ratio of the upper and 
lower bounds scales linearly with the aspect ratio k, rather than as k, 3 . In- 
deed, for conformal and anticonformal reflection-symmetric topologies (for 
which the wrapping numbers w acr about a given vertex have the same sign), 
one can show that the ratio is bounded by [L\ + L 2 y + L 2 Z )/L Z . It is not 
clear that for general prism topologies the k 3 dependence in Theorem 11.31 is 
optimal. 

An important question is whether within a given homotopy class the 
infimum Dirichlet energy is achieved. The homotopy classes 7iCp{h) are not 
weakly closed with respect to the Sobolev norm, so it is not automatically the 
case that the infimum is achieved. However, while a given 7iCp(h) may not be 



K_Ep{h) < inf E FO (n) < CK + K 3 E P (h) 



(1.15) 




(1.16) 
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weakly closed, it may still contain a local (smooth) minimiser of the Dirichlet 
energy Indeed, there is some numerical evidence and heuristics to suggest 
that, in the case of a rectangular prism, for the simplest topologies a smooth 
minimiser always exists, while for others a smooth local minimiser may or 
may not exist depending on the aspect ratios [TBI US], ft would be interesting 
to establish for which topologies there exist smooth local minimisers, also 
from the point of view of device applications. Such configurations would of 
course satisfy the bounds established here. 

There is an extensive literature on S^-valued harmonic maps with fixed 
(Dirichlet) boundary data; reviews are given in jHllS]- Problems related to the 
one considered here concern liquid crystal droplets [131122], in which one seeks 
configurations n on a three-dimensional region Q which minimise the elastic 
energy [Tl]. In case of tangent boundary conditions, there are necessarily 
singularities on the surface of Q, eg 'boojums' [21]; in a polyhedral domain, 
these singularities are pinned at the vertices. 

The remainder of the paper is organised as follows. Theorem ll.2l is proved 
in Section [21 and Theorem 11.31 in Section El modulo two lemmas concerning 
the explicit construction of and estimates for the representative prism con- 
figurations (Sections El and E]). In the Appendix it is shown that homotopy 
classes of continuous tangent unit-vector fields on P are classified by wrap- 
ping numbers. 

2 Lower bound for general polyhedra 

Proof of Theorem \l. 6 A In ^Tj\ we show that smooth n's are dense in 7iCp{h) 
with respect to the Sobolev W /< - 1,2 - l -norm. Therefore, it suffices to establish 
the lower bound ()1.12|) for n smooth. 

Let B e {y a ) denote the e-ball about v a , and let 



denote the domain obtained by excising these balls from P. Clearly 



Let x" denote the characteristic function of the sector S a C S 2 . It will 
be useful to introduce smooth approximations x a to such that x" has 
support in S a and satisfies < x" < X° ■ Then 



p e = p-u a (s e (v a )nP) 



(2.1) 




(2.2) 




(2.3) 
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Using the inequality jl] 

(Vn) 2 > 2||n*w||, (2.4) 

where n*u denotes the pullback of uj by n and || ■ || denotes the norm on 
forms induced by the standard metrics on R 3 and S 2 , we get that 



E(n)>2j2 I \\tf{ru)\\dV. 



(2.5) 



For each a, let £ CT denote a continuous piecewise-differentiable function 
on P with 

iidrn<i- (2.6) 

Then for arbitrary a, b, c, we have that 

||n*o;|| \dV(a,b,c)\ > (d£ a A nV)(a,b, c), (2.7) 
where dV is here regarded as the Euclidean volume form on M 3 . But 

dC A = d(C A n*(x a u)), (2.8) 

since d(n* (x° uj)) = n*d(x a uj) = (x CT co> is a two-form on S 12 ). Therefore, 

E(n)>2j2 ! ^(fAn*(f W )). (2.9) 

From Stokes' theorem, ()2.9j) implies that 

E(n)>2j2f rn*(f W ), (2.10) 

The boundary of P e consists of i) the faces of P with points in B € (v a ) removed, 
and ii) the intersections of the two-spheres dB € (v a ) with P. The latter, 
denoted by C e a = dB t (y a ) fl P, we call cleaved surfaces. Tangent boundary 
conditions imply that n*u; vanishes on the faces of P (since the values of n 
on a face are restricted to a great circle in S 2 ). Therefore, only the cleaved 
surfaces contribute to the integral in (|2.1()jl . We obtain 

E(n)>2j2J2 [ Wr"). (2-11) 
Jc a 

By the Bounded Convergence Theorem, we can replace x a by \ u in (j2.11j) . 
Taking the limit e — > 0, we obtain 

E(n)>2j]^r(v a )lim I n*( X ^). (2.12) 



From (jl.9|) . the integral over C e a yields A a times the wrapping number w acr , 
which depends only on the homotopy type of n. Thus, 

a a 

The remainder of the argument proceeds as in j^j. We note that (12.13)1 
holds for any choice of £ CT 's consistent with the constraints ()2.6j) . These 
constraints imply that 

|r(v a ) - T(v b )| < |v a - v 6 |. (2.14) 

Conversely, given any set of values £ a<J for which |£ a<T — £, bcr \ < |v a — v 6 |, we can 
find functions £ CT which satisfy the constraints ()2.6|) and assume these values 
at the vertices (for example, take £ CT (r) = max a (£ acr — |r — v a |)). Thus, we 
obtain a lower bound for Ep(h) in terms of the solutions of a finite number 
of linear optimisation problems, one for each sector, 

E'f(h) > A a I max ^ J^fV ] . (2.15) 

A simpler characterisation is provided by the dual formulation, 

EfHh) > 2 V A a [ min V |v a - v b \tt ab ' a ) , (2.16) 

c V a,b / 



where the Q ab,cr, s are constrained by 



Let us fix a. Without loss of generality, we can restrict VL ab,a to be non- 
negative and equal to zero unless w aa > and w ba < 0. Suppose first that 
the nonzero wrapping numbers are either +1 or — 1; by (jl.llj) there are an 
equal number, m say, of each. Therefore, the nonvanishing elements of VL ab > (T 
may be identified with an m x m matrix, which we denote by M. f|2.17|) 
implies that M is doubly stochastic. By a theorem of Birkhoff [2 , M can be 
expressed as a convex linear combination of permutation matrices. Then the 
minimum in ()2.16j) is necessarily achieved at an extremal point, ie for M a 
permutation matrix corresponding to a minimal connection. In this case, 



mm 

n ab >< 7 



K - v fc |ft a6 '< 7 = L(V a (h),M a (h)), (2.18) 
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where V c {h) (resp. N a (h)) contains vertices v a for which w aa equals +1 
(resp. —1). The case of general nonzero wrapping numbers values is treated 
by including v a with multiplicity \w aa \ in either V a (h) (for w aa > 0) or 
N a (h) (for w aa < 0). The same argument applies in every sector (there is a 
separate minimal connection for each a), and ()1.12j) follows. 

□ 

3 Upper bound for rectangular prisms 

Let P denote a rectangular prism centred at the origin of three-dimensional 
Euclidean space. We take the edges of P to be parallel to the coordinate 
axes and of lengths L x , L y , L z , oriented so that L x > L y > L z . It will be 
convenient to introduce the half-lengths 

lj = Lj/2 (3.1) 

(here and in what follows, the index j takes values x, y or z). Then the 
vertices of P are of the form 

v Q = (±l x ,±l y ,±l z ), (3.2) 

where the vertex label a designates the signs in ()3.2j) . 

Let O a C S 2 denote the spherical octant of directions about v a which 
are contained in P. Eg, for v a = (— l x , —l y , —l z ), O a is the positive octant 
{s e S 2 | Sj > 0}. The boundary of O a , dO a , contains the directions which 
lie in the faces at v a , and is composed of quarter-segments of the great circles 
about x, y and z. Let dO? denote the segment about j. 

Choose / so that < / < l z . Then v a + lO a is contained in P, so that 
v a + lO a is a cleaved surface. Given n e Cp(h), we can define a unit- vector 
field v a on O a by 

i/ a (s) = n(v a + /s). (3.3) 

Tangent boundary conditions imply that, for s G 80°: , u a (s) is orthogonal 
to j. Denote the set of u al s collectively by v. The wrapping numbers of n, 
and hence its homotopy type, are determined by v. i> is an example of an 
octant configuration, which we define generally as follows: 

Definition 3.1. An octant configuration i> with admissible topology h = 
{w aa } is a set of continuous piecewise- smooth maps u a : O a — > S* 2 satisfying 
tangent boundary conditions, 

sGSOf^ Ks)-j = 0, (3.4) 
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such that 

u*( X a uj) = w aa A a . (3.5) 



O" 

The Dirichlet energy of an octant configuration v on the octant O a is 
defined by 

Efa(") = / (Vi/°) 2 (s)dfi a . (3.6) 

J O a 

Here and in what follows, it will be convenient to regard Vi^(s) (the gradient 
of the jth component of u a ) as a vector in IR 3 which is tangent to O a at s. 
dQ a in ()3.6|) denotes the area element on O a (normalised so that O a has area 
7r/2). The Dirichlet energy on the octant edge dO°: is given by 

-I 2 ( d 



Here, s|(a) denotes the parameterisation of dOj by arclength (ie, angle) a. 
For example, if v a = (— l x , —l y , —l z ), 

s£(a) = cos ak + sin al, < a < vr/2, (3.8) 

where (j, k, I) denote a triple of distinct indices. 

By an extension of an octant configuration u, we mean a continuous, 
piecewise-smooth unit- vector field n on P such that n(v a + Is) = v a (s) for 
all s G O a . Obviously, if v has topology h, so has its extension n. We 
introduce the following notation: Given functions / and g on a domain W, 
we write / < g to mean there exists a constant C such that |/| < C\g\ on 
HA In this case, we say that / is dominated by g. 

Lemma 3.1. Let v be an octant configuration with admissible topology h. 
Then v can be extended to a continuous piecewise- differentiable configuration 
n G HCp(h) such that 

E(n) < k*L z ( Y,EU») + E ( E hM 1/2 J • ( 3 - 9 ) 

\ a aj / 

Theorem 11.31 is proved by constructing octant configurations whose Dirichlet 
energies on the octants and their edges scale appropriately with the wrapping 
numbers. These configurations are provided by the following: 
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Lemma 3.2. Given an admissible topology h = {w aa }, there exists an octant 
configuration v with topology h such that 

Y,Ei 2) {»)< x>n (3.10) 

a aa 
aj aa 

The proofs of Lemmas 13.11 and 13.21 which involve explicit constructions and 
estimates, are given in Sections |U and |S] respectively. 

Proof of Theorem MJA Given an admissible topology h = {w acr }, we choose 
an octant configuration i> with topology h as in Lemma 13.21 and extend it 
to a unit-vector field n on P as in Lemma 13.11 From the Cauchy-Schwartz 
inequality (|3.1H) implies that 

E%w) 1/2 ^Ei wa<T i- ( 3 - 12 ) 

aj aa 

Together, (|3.9jl . (|3.10|) and ()3.12|) provide an estimate for E(n), and therefore 
an upper bound for Ep { (h), 

Ef(h) < E(n) < k 3 L z ^ \w aa \- (3-13) 

aa 

From Theorem 11.21 a lower bound for Ep (h) is given by 

E P (h) = J2^HP a (h),N^h)) (3.14) 

a 

(A a = 7r/2 for a rectangular prism). The minimum distance between vertices 
of P is L z . As the number of elements olP a (h) (and of N a {h)) is | Y2 a l^ a °1) 
it follows that 

L(P a (h), N a (h)) > \L Z \w aa \- (3.15) 

a 

From (|3.13|) and (|3.15jl . we conclude that 

E^(h) <n 3 E P (h). (3.16) 

□ 

We remark that the octant configurations of Lemma 13.21 must be chosen 
with some care, as the following example illustrates (details may be found 
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in For simplicity, take P to be the unit cube. In the (positive) octant 

about v a = (— |, — |, — |), let 

v a (9,(f>) = (sin a cos /3, sin a sin /3, cos a), (3-17) 

where < 0, <\> < ir/2 and a = (AM + 1)0, /3 = (AN + 1)0 for integers M 
and AT. Given (x, y, z) & P with x,y,z < 0, let (0, 0) denote the polar angles 
of (x,y,z) with respect to v a , and let n(x, y, z) = u a (9,(j)). We define n 
elsewhere via n(±x, y, z) = n(x, ±y, z) = n(x, y, ±z) = n(x, y, z) (so that n 
is a reflection-symmetric configuration [TT)1 118j). Then n is continuous and 
satisfies tangent boundary conditions. Denote its homotopy class by Hmn- 
It is straightforward to compute the wrapping numbers (it turns out that 
they scale linearly with M and N) , and, from Theorem II. 2[ to obtain the 
following lower bound: 

E P (h MN ) = (2 max(M + 2N, 2M + N) (3.18) 

It turns out that E(n) can be evaluated exactly as a finite sum of Appell 
hypergeometric functions. For large M and N, the energy is given asymp- 
totically by 

E(n) ~ 4v / 3((4M + l) 2 + \ In M(AN + l) 2 -. (3.19) 

Clearly E(n) does not scale linearly with M and N, so is not dominated by 
the lower bound Ep(hMN)- 

4 Extending octant configurations 

Let us specify the geometry of the prism P in greater detail. Let 

m^=v*-^j (4.1) 

denote the midpoint of the edge through v a along j (here, is the jth com- 
ponent of v a ). Let C a denote the triangle whose vertices are the midpoints 
of the edges coincident at v a , 

C a = {r = + r y m a {y) + T z m? z) | Tj > 0, ^ Tj = 1}. (4.2) 

3 

We call C a a cleaved face (see Figure [TJ). c E C a satisfies 

C a ■ c = 2, (4.3) 
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(b) 



(c) 





(d) 



Figure 1: (a) The cleaved plane C a (b) The pyramid X a (c) The pyramid 
Y a (d) The truncated face F^ T (e) The pyramid Z^ T 



where 

C a = H)- 1 ,^)- 1 ,^)- 1 ). (4.4) 

is an (unnormalised) outward normal on C a . Let h denote the distance from 
C a to the origin. Then 



-L < h 



< 2L 



(4.5) 



Let F^ T , where r = ±1, denote face of the prism which lies in the plane 
{ r j = r (j'}- Let F^ T C F^ T denote the rhombus whose vertices lie at the 
midpoints of the edges of F jV . We call F^ T a truncated face (see Figure HJ). 

We partition P into three sets of pyramids, denoted X a , Y a and Z JT . X a 
and Y a have the cleaved face C a as their (shared) base. X a has its apex at 
v a , while Y a has its apex at the origin. Z^ T has the truncated face F^ T as 
its base and its apex at the origin. Every point of P belongs either to the 
interior of just one of these pyramids or else to the boundary between two 
or more of them (see Figure . 

In the proof of Lemma I3.1[ we define n, the extension of the octant 
configuration is, to be constant along rays from C a to the origin (apart from 
a small neighbourhood thereof) and along rays to v a . We then show that 
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the energy of n in X a and Y a is proportional to E%,(y). The construction 

of n in Z^ T is more complicated. On dF^ (the boundary of the base), n 
is determined by v, and in the interior of F^ T , we define n by a simple 
interpolation which respects tangent boundary conditions. In the interior of 
Z^ T , we do not take n to be constant along rays from the apex to the base, 
as this would give rise to an energy proportional to ^2 v a eF jr E^j(u), which, 
for the octant configurations of Lemma 13. 2\ would scale as the square of the 
wrapping numbers. Instead, along such rays, and over a distance 

°= (* E > ( 4 - 6 ) 

V v a eF^ T / 

n is rotated toward the normal j. This leads to an energy in Z jV proportional 
to l/cr. 

Proof of Lemma \H. 11 Given an octant configuration u, we define an exten- 
sion n on the pyramids X a , Y a and Z^ T (Steps 1-3) with Dirichlet energies 
-Ex^n), E ya (n) and E Z jr{vL) bounded as follows: 

E X a(n) < ±kL z E? 2) (v), (4.7a) 
E Ya (n) <K 3 L z Efo(v), (4.7b) 

E Z ^)<^L Z E (^(i) J ») 1/2 , (4.7c) 

Then 

E(n) =E(^(n)+ J By»(n)) + E^(n) 

a jr 

< k 3 L z [ E^( 2 )(^) + E (£(D>)) ) • ( 4 - 8 ) 



To ensure continuity (Step 4) we modify the construction of n near the origin 
while preserving the bounds ()4.7I) . 

Step 1. Construction in X a From (|4.2j) and ()4.3j) . points in X a are of the 
form v a + rs, where s G O a and < r < r a (s) with 

r a (s) = (4.9) 
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r a (s) is the distance from v a to C a along s. The maximal distance is half 
the length of the longest edge, so that 

r a (s) < r a (x) = l x . (4.10) 

We define n in X a by 

n(v a + rs) = j/*(s), s G O a , < r < r a (s). (4.11) 
Then, from O and (ICTIj) . 

E xa (n):= [ {VrvfdV= [ r a (s)(W a (s)) 2 d{l a 
Jx a Jo a 

<l x El 2) (u)<\KL z E a {2) (v), (4.12) 

as in ()4.7ajh 

Step 2. Construction in Y a Points in Y a (excluding the origin) are of the 
form Ac, where c G C a and < A < 1. We define n in Y a by 

n(Ac) = n(c). (4.13) 

Note that n(c) is fixed in Step 1, since C a belongs to X a as well as Y a . 

To estimate .Eya(n), we resolve Vn into components tangent and normal 
to the cleaved face C a , 

Vn = V t n + V„n, (4.14) 

so that (C a ■ Vf)n = and (t • V n )n = for t tangent to C a . From ifPTty . 
Vfn(Ac) = A _1 V(n(c). Therefore, 

(V,n(Ac)) 2 = ^(V,n(c)) 2 . (4.15) 

To estimate (V n n) 2 , we note that (c • V)n(Ac) = (n is constant along rays 
in Y a through the origin) and resolve c into components c t and c n tangent 
and normal to C a to obtain 

(V n n(Ac)) 2 < |^(V,n(Ac)) 2 < ^^(V t n(c)) 2 . (4.16) 
Since (V 4 n(c)) 2 < (Vn(c)) 2 , (jUTjj) and (jj~IE}) together give 

(Vn(Ac)) 2 = ^^(Vn(c)) 2 . (4.17) 
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Clearly |c| < l x while |c n | is just h, the distance from C a to the origin, and 
h > 2/ z /3 (cf dHU)), so that 

(Vn(Ac)) 2 < ^ 2 l(Vn(c)) 2 . (4.18) 

The volume element on Y a is given by 

dV = h\ 2 d 2 cd\, (4.19) 
where d 2 c is the Euclidean area element on C a . Since h < 21 z , 

E Ya (n)= [ (Vn) 2 dV < \k 2 I z [ (Vn) 2 (c) d 2 c. (4.20) 

JY a 2 J C a 

Letting s = (c — v)/|c — v|, we can write the preceding as an integral over 
O a . We have that 

d 2 c = |s ^| 2 /|ca| ^ a , (Vn(c)) 2 = ^(W«(s)) 2 , (4.21) 

and 

|C a | 3// 2 

so that (|4.2()jl becomes 

27 f 

E Ya (n) < ——k?I z / (Vu a (s)) 2 dQ a < k 3 L z E%(v), (4.23) 
^ Jo a 

as in lj£7Kjl . 

S'tep 5. Construction in Z JT To simplify the discussion and the notation, 
let us fix our attention on the top face of the prism, with j = z and r = 1, and 
henceforth drop the designation jr, writing Z for Z^ T , F for F^ T , etc, in what 
follows (the other faces are handled similarly). OF may be parameterised as 
R(0) = (R((p) cos(p,R(4>) sin0, l z ), where 

R(d>) = J J « , , - (4.24) 

Ly\ cos 0| + l x \ sm0| 

and is the polar angle about the z-axis. On dF (which is also contained in 
the cleaved faces), n is defined in Step 1. It follows that n is continuous on dF 
(including the midpoints of the edges of F, which belong to two cleaved faces, 
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as v has an admissible topology) and satisfies tangent boundary conditions 
there. Tangent boundary conditions imply that 

n(R(0)) = e cos(6(0))y + sin(6(0))x, (4.25) 

where 0(0) may be taken to be continuous and piecewise smooth, with 
6(27r) — 6(0) equal to a multiple of 2n. Since v> has an admissible topology 
we must have ®{2n) = 0(0). e = ±1 can be chosen so that 

6(0) = 6(27r) = 0. (4.26) 

We introduce polygonal cylindrical coordinates (h, £, 0) on Z, defined by 

x = h^R{(p) cos0, y = h£R(<f>)sm<f>, z = l z h, 

0<h<l, < f < 1, < < 2tt. (4.27) 

£, the radial coordinate, is scaled to equal 1 on the sides of Z and along 
the z-axis. Then let 

N (2) (£,0) =ecos(£e)y + sin(£9)x. (4.28) 

N(2) gives a continuous extension of n to the interior of F which satisfies 
tangent boundary conditions. A continuous extension to the interior of Z is 
given by 

N(/i, f , 0) = cos 7N (2 ) + sin 7Z, (4.29) 
where 7 = j(h, is given by 

^..(^.(izi)- ^.{j. : <; ; (4 ,o, 

and < cr < 1. Thus, 7 vanishes on the boundary of Z and has a constant 
value, 7r/2, at interior points sufficiently far from the boundary, a, which 
determines how far, will be specified below. We define n on Z as 

n(x(h, e, 0), e, 0), e, 0)) = n(/i, e, 0). (4.31) 

It is readily checked that ()4.31|) agrees with (j4.13|) at points on the boundary 
of Z except at the origin. 

The energy of n in Z is given by Ez{p) as follows. In terms of the 
coordinates (h, £,0), we have that 

E z (n) = [ {VnfdV 
J z 

= dh dU #|(Ve x V0) • V/ip^NftV/i + N^Ve + IVWO 2 , 
Jo Jo Jo 

(4.32) 
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where N h = dN/dh, N c = <9N/<9£, etc. From (l4~2Kl) and we get that 

N/, = 7/,,(cos7Z - siri7N (2 )), 

= 7^(cos7z — siri7N( 2 )) + e0cos7N( 2 ) x z, 

= e£0' cos 7N (2) x z. (4.33) 

From 1)4.27)1 . we get that 

vh= (o,o,i/y, 

V£ = £ ( (i2 cos + i2' sin 0) / (pR) ,(Rsm<f)-R' cos 0) / (pR) ,-l/(l z h)), 
V0 = (- sin 0, cos 0, 0) /p, (4.34) 

where p = (x 2 + y 2 ) 1 ^ 2 = h£R. Straightforward calculation then gives an 
expression for Ez(ji) of the form 

E z (n) = y^St = V f dh f di r d4I u (4-35) 
i=i »=i ^0 Jo Jo 

where the integrands for the separate contributions £j are given by 

h = h COS 2 7^0' 2 , h = ^ 7^, h = -2y lhlii eR\ 

TP 

h = -2e/ z cos 2 7^— 00', 
R 

h = (u ( 1 + jP) + X R2 ) ( C ° S2 + 7 ^ • (436) 

We consider these contributions in turn. 

Concerning S\, since cos 2 7 vanishes for < £, h < 1 — a (cf (14. 30)1 ). it 
follows that 

1-2-K 

£ x < 2l z a I ,2 (0) dcj). (4.37) 
Jo 

The integral of /2 (0) can be related to the Dirichlet edge energies E^ z {y) 
for the vertices v a which lie on F. For convenience, label these anticlockwise 
by a = 0, 1, 2, 3 so that 

n(R(0)) = i/ a (s^(a!(0)), (a - l)vr/2 < < avr/2, (4.38) 

where s"(a) parameterises the octant edge 90" as in ()3.8)1 . and a((f>) is given 
by 

{{IJlyf tan 0, O<0<7r/2or7r<0<37r/2, 
tana = < „ (4.39) 

|(^//x) 2 tan0, vr/2 < < tt or 3tt/2 < < 2tt. 
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Figure 2: The angles a and 



(4> is the angle with respect to the centre of F and a, with < a < tt/2, the 
angle with respect consecutive vertices. ()4.39j) gives the elementary relation 
between them. See also Figure El) Recalling ()4.25|) . we get that 



d 



d 



da 



a=a(tf>) 



f da\ 



It follows that 

-2tt 







^W=ErO«»] [fa ] da - 



v a £F 

From 1)4.39)1 one has that 



da 



I 2 

< JL < K 2 
- 72 - K 



Therefore, 



re' 2 w#<« 2 5]i? ( » 1)2 (,) 

JO v a eF 

Substituting into (J4.37)) . we get 



Next we consider 



Si 



v a eF 



1 r 2 pi /-27T 



I) ^2 JO 



(4.40) 



(4.4i; 



(4.42) 



(4.43) 



(4.44) 



(4.45) 
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From ()4.24j) . R(<p) < l x , and from (|4.30jl . |7/J vanishes for h < 1 — a and is 
bounded by n/{2a) for h > 1 — a. Therefore, 



to < — > K U 



4 La 



a 



(4.46) 



We estimate £3 similarly; noting that |7g| vanishes for £ < 1 — a and is 
bounded by ir/ (2a) for £ > 1 — a, we obtain 



£ 4 is given by 



a 



(4.47) 



£4 = —2el z / dh I d£ £ cos 7 
'0 </o 



2tt 



(4.48) 



We consider the ^-integral first. From the Cauchy-Schwartz inequality, 

1/2 







2n d/ 

06'^ dcj) 
R 



< 



■2n 



1/2 / r 2n / p/\ 2 



(4.49) 



In the first factor on the right-hand side, note that, from ()4.25|) and (j4.38|) . 



,tt/2 

10(0)1 < E / 



^«(a)) 



da, for < </> < 2tt. 



(4.50) 



The Cauchy-Schwartz inequality then implies that 

6 2 (0) < 2tt E forO<0<27r. 



(4.51) 



In the second factor on the right-hand side of (|4.49|) . we have that \R'/R\ < k 
(cf flPg) ), so it follows from (IQ3l that 







2?T / T)l 

9 ' 2 fi 



< 



* 4 E 



(4.52) 



v«eF 



Substituting (|43TJ) and (l4~H2l into (IP9l . we get that 



27T JV 



(4.53) 



21 



We substitute ()4.53|) into ()4.48|) and recall that cos 2 7 vanishes for < £, h < 
1 — a to get that 

i54i<«vE^( a i)»- ( 4 - 54 ) 

Finally, we consider 

£ 5 = j x d/» jf 1 de jT rf0 (cos 2 7 e 2 + 7 |) ( x + + f • 

(4.55) 

Let us first estimate the ^-dependent terms. From (|4.24|) . we have that R < 
kI z and \R'/R\ < k, while f|4.51|) provides a bound for 6 2 (0). Substituting 
into (|4.55|) . we get that 

£ 5 < 2nl z [ dh [ dd2n cos 2 7 E + 7 | J (£(1 + k 2 ) + £V)) . 

(4.56) 

Recalling that cos 2 7 and 7? vanish for < £, /i < 1 — a while 7| is bounded 
by 7r/(2cr), we obtain the bound 



S 5 < Airl z a ( 2tt E £(\ )z (i/) + ^ ) (1 + 2« : 

v a £F 



2\ 



<^fe%M H 2 lA (4.57) 



We substitute the estimates (jQ7j) . (jOfijl . (jQ7|l . (l4~olll and (PTKoT) for ^ 
into (|4.35J) . and take 



-1/2 

(4.5? 



We can verify that a < 1 as follows: Tangent boundary conditions imply that, 
for a multiple of 7r/2 (ie, for R(0) belonging to an edge of F), ©(</>) = 7r/2 
mod 27T. Then (14. 51 J) implies that X] v «eF s nb > 7r/ 8 > so tnat a < y/8/n < 1. 
Then (ET37I) and (ET3K1) give 



\v a GF / v a eF 



1/2 , (4.59) 



as in (|4.7cjl . The same estimate may be carried out for the other pyramids 
Z^ T . Different lengths lj appear as appropriate, but since ratios of lengths 
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are bounded by k, f!4.59|) holds generally. 



Step 4- Continuity As defined, n is continuous everywhere except at the 
origin. Here we modify n in a small neighbourhood of the origin to remove 
the discontinuity while preserving the estimate (|4.7j) . 

From the definitions ()4.13|) and ()4.31|) of n in Y a and Z^ T , n(r) is radially 
constant in the ball B$ about the origin of radius 5 = (1 —a)l z . Let 7 : S 2 — > 
S 2 denote the restriction of n to dB$, given by 

7(s) = n(5s), s G B s . (4.60) 

Let 

E ( 2)h)= [ (v 7 ) 2 ^ (4.6i) 

Js 2 

denote the Dirichlet energy of 7. Then, for < e < 1, the energy of n in the 
e<5-ball about the origin is given by 

E B Jn) = e6E {2) ( 7 ). (4.62) 

7 is piecewise smooth, and, since v has an admissible topology, of degree 
zero. It follows that 7 is smoothly homotopic to a constant map. Let Y s 
denote a homotopy, so that Ti(s) = 7(s) and r (s) = s . Let g be the 
unit-vector field on Bs given by 

g(rs) = T r/S (s), (4.63) 

and let E Bg (g) denote its Dirichlet energy. For < e < 1, we define g e to be 
the unit-vector field on B e $ given by 

g e (r) = g(r/e). (4.64) 

Then the energy of g e in the e<5-ball about the origin is given by 

E B Jg e )=eE Bs (g). (4.65) 

For any < e < 1 we can redefine n on B e s, taking it to be g e there, and 
leaving n unchanged elsewhere. It is clear that n as redefined is continuous 
and piecewise smooth. From ()4.62)1 and ()4.65j) . the redefinition changes its 
energy by e(EB s (g) — 5E(2)(l))- As e can be made arbitrarily small, the 
estimates ()4.7|) remain valid. 

□ 
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5 Constructing octant configurations 



For definiteness, we consider the configuration on a particular octant, namely 
the positive octant about the vertex (—l x ,—l y ,—l y ); the treatment for the 
other octants is analogous. To simplify the notation, we will drop the vertex 
label a. Hence, throughout this section, we write O = {s £ S 2 \sj > 0} for 
the positive octant (instead of O a ), u : O —>■ S 2 for the configuration on O 
(instead of v a ), w a instead of w aa , etc. With reference to (jH.fi j) and (j3.7j) . we 
let 

E<${y)= [ (W) 2 (s)dft, (5.1) 
Jo 

r /2 ( dv\ 2 

%W = J [^) Ma)) da, (5.2) 

where <El denotes the area element on O and 

Sj(a) = cosak + sinal (5.3) 
denotes the parameterisation of dOj. Lemma f3.2l follows from showing that 

Ew{v)<52\w% (5.4) 

a 

Y,E {1)j (u)<J2\™ a \ 2 , (5-5) 

as analogous relations hold for the other octants. Before establishing (j5.4)l 
and (j5.5j) in Section 15.31 we first review the topological characterisation of 
octant configurations fSection l5.1j) and their representation by complex func- 
tions, particularly conformal representatives ( Section l5.2j) . 

5.1 Topological characterisation 

As discussed in [21] (in the context of general convex polyhedra) and in 
[Td^ ITS] (for a rectangular prism), the homotopy class of u : O —>■ S 2 may 
be characterised by certain invariants, namely the edge signs, denoted e = 
(e x , e y , e z ), kink numbers, denoted k = {k x , k y , k z ) and trapped area, denoted 
Q. Here we recall the definitions of these invariants and some relevant results 
for prisms; details may be found in the references. 

Tangent boundary conditions imply that v(j) is parallel to j; the edge 
sign ej determines their relative sign, ie 

= ejl (5.6) 
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Tangent boundary conditions also imply that along Sj(a), v takes values in 
the (fcZ)-plane. The integer-valued kink number kj counts the number of 
windings of u in this plane relative to the minimum possible winding (a net 
rotation of ±7r/2), for which kj = 0. The trapped area Q is the oriented area 
of the image of u, given by 

n = [ u*cu. (5.7) 
Jo 

For a rectangular prism, the sectors are octants of S 2 labeled by a triple 
of signs o = (a x , a y , cr z ), with S a = {s G S 2 | sgn (s • j) = <jj}. The wrapping 
numbers w a may be expressed in terms of (e, k, Q), as follows |2~T) I18j: 

3 

where S ae equals one if a = e and is zero otherwise. Note that ()5.8|) implies 
that 

n = —2n'^^ajkj — e*7r/2 mod An, (5-9) 

3 

where 

e* = e x e y e z . (5.10) 
(|5.8j) may be inverted to obtain (e, fc, f2) in terms of the wrapping numbers, 

e t = — a r a s w a , for r, s, t distinct, (5-11) 
k r = \ (x r w CT + |e r e*, (5.12) 

fi = Ei wff - ( 5 - 13 ) 

(Similar relations are described for a general polyhedron in the Appendix.) 

A topology for a prism configuration n may be specified as a set of edge 
signs, kink numbers and trapped areas for each vertex. The conditions for 
the topology to be admissible (ie, realisable by a configuration continuous 
away from the vertices) are readily expressed in terms of (e,k,Q); pairs of 
edge signs associated with a single edge must be compatible, the absence of 
surface singularities implies sum rules for the kink numbers on each face of 
the prism, and the absence of interior singularities implies a sum rule for the 
trapped areas. 

We say that an octant topology (e, k, O) is conformal (resp. anticonformal) 
if every nonzero wrapping number is negative (resp. positive). From (|5.8jl . 
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one can show that (e, k, fl) is conformal (resp. anticonformal) if and only if 
Q < —fl_(e,k) (resp. Q > Q + (e,k)), where 



^(6,^=2^^1^1+2 




if e^Cjkj < for all j, 
otherwise. 

tu.fet) = ar^w - 2 - 1 T r if ,!* eA > for a " h < 5 - 14 > 

L — ' 1 — - otherwise. 

If (Q, e, k) is an octant topology, then Q differs from Q + (e, k) or fi_(e, fc) by 
some multiple of 4n (ie, some number of whole coverings of the sphere). If 
Q = fi±(e, k), then (e, fc, f2) has at least one wrapping number equal to zero. 



5.2 Conformal configurations 

v : O —>■ S 2 can be represented by a complex function F(w, w) in the stan- 
dard way via the stereographic projection S 2 C* (C* is the extended 
complex plane), 



^ x I ^x ItS 



1 + s z 1 + s z 



v x + iv, 



s . 



(5.15) 



The domain of F is the positive quarter-unit disk (the image of O under the 
projection), 



Q = {w e C| < Rew < 1, < Imw < 1, |w| < 1}. 



(5.16) 



Letting Wj(a) denote the projections of the parameterised boundaries Sj(a) 

(5.17) 



of dOj, we have that 



w x (a) = tan(a/2) } w y (a) = i(l — tan(a/2)), w z = e 
A standard calculation gives 



(W) 2 = 4 



I- 1 w 


2 + l 


F 

1 w 


2 


(1 + 


F 


2 ) 2 



(5.18) 



so that (cf (jnHJ) and (Q) 



£(2)(") 





2 + |^F| 2 


(1 + 


F 


2y 



tt/2 



d w, 



\d^F\ 2 + \d w F\ 
(1 + |F| 2 ) 2 



{wJa)) 



dvjj 
da 



Oi, 



da. 



(5.19) 
(5.20) 
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Of particular interest are configurations for which F is conformal (a func- 
tion of w only) or anticonformal. For definiteness we will consider the con- 
formal case, and write F(w,w) = f(w). If / has a meromorphic extension to 
the extended complex plane, then tangent boundary conditions imply that 
/ is real when w is real, / is imaginary when w is imaginary, and |/| = 1 
when \w\ = 1. It follows that if w* is a zero of /, then so are -w t and ±tu*, 
while ±l/iw* and ±l/iu* are poles. Therefore, / is in fact a rational function 
of the form 

/ = Xw n ABC. (5.21) 

Here, A contains the real zeros and poles of /, B the imaginary zeros and 
poles, and C the strictly complex zeros and poles, n is an odd integer giving 
the order of the zero or pole at the origin, and A = ±1 is an overall sign. A, 
B and C may be written explicitly as 

A{w) 

B(w) 
C(w) 

Here, a is the number of real zeros and poles in Q (excluding the origin), b 
the number of imaginary zeros and poles in Q (excluding the origin), and c 
the number of strictly complex zeros and poles in Q. rj, with < Tj < 1, 
denotes the real zeros (pj = 1) and poles (pj = —1); similarly, isk, with 
< Sfc < 1 and = ±1, denote the imaginary zeros and poles, and ti, with 
< \ti\ < 1, Ret/, Imti > and t\ = ±1 denote the strictly complex zeros 
and poles. 

In terms of these parameters, the edge signs, kink numbers and trapped 
area of conformal configurations are given by ^Hl EB1 

e« = A(-l)°, e, = Al-l)'!-!)!"- 11 ' 2 , e 2 = sgn n, (5.23) 



11 I r 2 w 2_ 1 J 

nf ^- g) 
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\fe=i 

fc, = -\{-l) a e x \J2i-lYpj + |(1 - (-l) a )e^ , (5.24) 

a be 

k z = \ (e x e y - n) - \^Pj - |^ afc ~ 

j=i k=\ i=i 

Q, = -\{\n\ + 2(a + b) + 4c)tt. (5.25) 

The expressions for the edge signs follow from evaluating / at 1, i and 0, 
while the formula for Q follows from noting that 8Q is just (— 47r times) the 
degree of /, which for a meromorphic function is the number of its zeros 
counted with multiplicity. The formulas for the kink numbers require a bit 
more calculation; details are given in [TH] . 

It is easily checked that a conformal configuration f(w) has a conformal 
octant topology (/ is orientation-preserving, which implies that the nonzero 
wrapping numbers are negative). Therefore, octant topologies which are 
neither conformal nor anticonformal, ie (e, k, Q) for which — Q_(e, k) < Q < 
—Q + (e, k), cannot be realised by F = f(w) or F = f{w). In [THj we establish 
a converse result, namely that evey conformal (resp. anticonformal) octant 
topology has a conformal (resp. anticonformal) representative. 



5.3 Proof of Lemma 13.21 

Proof. Given an admissible topology h for P, let (e, k, Q) denote the associ- 
ated octant topology on O with wrapping numbers w a '. We construct v, or 
rather its complex representative F(w,w), with topology (e, k, ft) in Step 1. 
We establish the estimate (|5.4|) for E( 2 ){v) in Step 2. We then show that 

£(D»< 1 + (5-26) 

for j = z (Step 3) and j = x,y (Step 4). From (|5.9jl . Q ^ 0. (j5.13j) then 
implies that the wrapping numbers cannot all vanish, so that, from (|5.12|) . 

£l + N 2 <]>>1 2 - (5.27) 
The bound Q for V . E {1)j {v) then follows from flOgj) and (I5~2TI) . 



Step 1. Definition of F In general, the octant topology (e, k, fl) is neither 
conformal nor anticonformal. We will take F to be conformal outside a small 
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disk in Q with conformal topology (e,k,Q*). Inside the disk, F is made to 
cover the complex plane (Q — Q*)/(47r) times; this will ensure that F has the 
required topology. 
Let 



fi* = -2tt 



^l^l + l-fe^ . (5.28) 



Using (|5.14jl , one can check that f2* is equal to either — fi_ (e, k) or — (e, fc) — 
47r; in either case, (e, k, fi*) is a conformal octant topology. Denoting its 
wrapping numbers by w°, it follows that wl < and, from (j5.8J) and ()5.14|) . 
that at least one of its wrapping numbers either vanishes or is equal to —1; 

\wl°\ < 1 for some a , (5.29) 

An explicit conformal representative f(w) with octant topology (e, fc, f2*) 
of the form (|5.21|) and ()5.22|) is obtained by taking 

n = (2 - e*)e 2 , A = e^, 
a = 2|/c^|, pj = -(-l) j e x sgnk y , 1 < j < a, 
b = 2\k x \, a k = -(-l) k e y sgnk x , 1 < k < b, 

i / n I — sgn k z , I < \k z \, 

c=|* a | + i(l-e.), n=\ -' ' . (5.30) 

as can be verified from ()5.23|) - (J5.25|) . (The reason for introducing fi* - 
we could use — f2_(e, k) instead - - is that conformal representatives for 
(e, k, — Q-(e, k)) entail several special cases.) Let 

Aixm = n-n,(e,k). (5.31) 

Then, from ()5.9j) . m is an integer, and from (J5.8)) . 

W CT = w ° + m . (5.32) 

If m = 0, we let 

F(w,w) = f(w). (5.33) 

Otherwise, let Wo denote a regular point of / in the interior of Q, and let 
D € (w ) denote the open e-disk about w . Choose e sufficiently small so that 
D 2e (w Q ) is contained in Q and contains no poles of /. Let 

s{w,w) = (5.34) 
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(so that s varies between and 1 as \w — Wq\ varies between e and 2e). Then 
for m > we define 

{f(w), \w - w \ > 2e, 

sf(w) + (1 - s)(f(w ) + (w- w ) m ), e < \w - w Q \ < 2e, 
/(ty ) + e 2m (u? - w y m , \w-w \<e, 

(5.35a) 

while for m < we define 

{/(w), |tu - iu | > 2e, 

s/(u>) + (1 - s)(f(w ) + (w- w )-' m ), e < \w - w Q \ < 2e, 
f(w ) + e~ 2m (w — w ) m , \w — w \ < e. 

(5.35b) 

F coincides with / on dQ, so that F has the same edge signs and kink 
numbers as /, namely e and k. Let us verify that F has trapped area Q. For 
m = this is automatic. Otherwise, for definiteness, suppose that m > 
(the case m < is treated similarly). From (|5.7|) and ()5.15|) one can show 
that 

Q(F)= [ 4 |( ^f '"J^f |2 rfV (5.36) 



Divide the domain of integration as in ()5.35j) . The contribution from |ty — 
iOo| > 2e is, to 0(e 2 ), just the trapped area of /, namely fl*(e, fc). Con- 
sider next the contribution from D e (wo). From (|5.35|h F(D t (wo)) covers the 
extended complex plane m times with positive orientation, apart from an 
e m -disk about f(wo). It follows that its contribution to the integral in ()5.36|) 
is, to within 0(e 2m ) corrections, 4nm. The remaining contribution, from the 
annulus -D 2e (wo) — D e (wo), is 0(e 2 ). This is because the area of the annulus 
is 0(e 2 ), while the integrand in ()5.36|) may be bounded independently of e 
(since, by assumption, / has no poles in _D 2 e(wo))- Since the trapped area is 
an odd multiple of 7r/2 (cf (|5.9|l ). it follows that, for small enough e, F has 
trapped area f2*(e, k) + 4irm = Q. 

Clearly the topology of F does not depend on the positions of the zeros 
and poles of /. As will be evident in Step 2 below, neither does the octant 
energy E( 2 )(v), at least to leading order in e. However, the edge energies, 
E(i)j(v), do depend on the positions of the zeros and poles. As will be evident 
in Steps 3 and 4, to obtain good control of the edge energies, the a real and b 
imaginary zeros and poles in Q should be kept away from the origin, the unit 
circle and each other, while the c strictly complex imaginary zeros and poles 
should be kept close to the unit circle and away from the real and imaginary 
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axes and each other. Anticipating these requirements, in (|5.22j) we take 



1 j lk ( 1 \ 1/2 . a, 1 I 



J 4 ' 2a' 4 ' 26' ^ c+lj " ' vr 8 ' 4(c+l)' 

(5.37) 

These imply the properties 

1/4 < rj < 3/4, r i+1 - r,- = l/(2a), (5.38a) 
1/4 < s k < 3/4, s k+1 - s k = 1/(26), (5.38b) 
tt/8 <ai< 3tt/8, 1 - |t/| 2 = l/(c + 1) (5.38c) 

which will be useful in what follows. We note that, with rj and pj as given 
in ()5.3()j) and ()5.37j) . the real zeros and poles of / alternate along the interval 
(0,1]; similarly, with s k and a k as given in ()5.3()j) and ()5.37j) . the imaginary 
zeros and poles alternate along (0, i]. 

Step 2. Estimate of E( 2 )(v) The expression (|5.19jl for E^iy) differs from 
the expression (|5.3fij) for Q only in the relative sign of the w- and w)- derivative 
terms. Arguing as for ()5.36j) . we see that, to order e 2 , the contributions to 
E(2)(v) from Q — D 2e (w ), D 2e (w ) — D e (w ), and D e (w ) are, respectively, 
and 47r|m|. Therefore, 

E (2) (u) < + 47r|m| = £ |<|| + 47r|m|, (5.39) 

cr 

where we have used (I5.13j) and the fact that w° < 0. From ()5.32j) . \w°\ < 
\w a \ + \m\. Also, from (El and (fQ2|) . \m\ = \w a ° - w%°\ < \w ao \ + 1 < 
2 J2tr Substituting these results into ()5.39jh we get that 

E (2) (v) <J2\ wa \- ( 5 - 4 °) 

verifying 1)5. 4|) . 

Step 5. Estimate of E [l)z (v) On <90 z , F = f and |/| = 1. Also, from (jQ7jl . 
Idu^/dal = 1. Therefore, from (|5.20|) . 



Below we show that 



/■t/2 




/o 


/ 



2 



(e iQ ) da. (5.41) 



y (e fa )<l + |U (5.42) 
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which then yields the required estimate ()5.26|) for j = z. 
To verify ()5.42|) . we note that, from ()5.2H) . 



n 



A' 
~A 



B' 



B 



C 



c 



(5.43) 



From the expression ()5.22aj) for A and the positions ()5.38a)) of its zeros and 
poles one calculates that 



A , a / 2 / 



r 2 j-i + r 2 j 



+ 



r 2 j-i + r 2J 



(w 2 - rlj_ x )(w 2 - r 2 j) {{rjj^w 2 - l){r 2 2J w 2 - 1), 

(5.44) 

From ()5.38aj) . |r 2 j + r 2 j_i| < 3/2 while, for w G dO z , we have that \w 2 — 
r 2 |, \r 2 w 2 — 1| > 1/16. Therefore, 



.4' 



.4 



a/2 

<5><1, wedO z 

< * n 



(5.45) 



j=i 



A similar calculation (cf ()5.22b|l and ()5.38bj0 shows that 

B' 



B 



(w) < 1, we dO z . 



(5.46) 



It remains to estimate C / C. Without loss of generality, we may assume 
that c> (otherwise, C = 1 and C = 0). From (j5.22cjl and (j5.38cjl . 



a 



2e 



■— ± 



(i + N 2 ) + (i + l^l 2 



\ e 2ia _ £212 



e 2ia _ p|2 



(5.47) 



Write the denominators (j5.47j) as 



2ia t2\2 



tf\ 2 = (l-\t l \ 2 y + 4\t l \ 2 sm 2 (a l + a) 

£ 2ia _ t 2|2 = (1 _ | tj | 2) 2 + 4 | ti |2 ^2^ _ q) 



(5.48) 



From (j5.38cjl . one has that 1 - \U\ 2 = l/(c+ 1) while, for < a < tt/2, we 
have that n/8 < a + on < 7n/S and — n/2 < a — a t < tt/2, so that 



le^-tfl 2 > sin 2 (7r/8) 

2ia ±2\2 - ' 



ttr > 



; c + i) 2 7T 2 



+ —^{oii - ay 



(5.49) 
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(we have used (7r/2) sin \x\ > \x\ for \x\ < vr/2). Substituting f)5.49|) into 
(fOTj) and using a, = vr/8 + ln/(A(c + 1)) (cf lET37t ) and 1 + |i,| 2 < 2 (cf 
()5.38c|l ). we get that 



C" 



^E( — *+V + 



^ Vl + |(^-(c+l)(4a/7r-i)) 2 sm 2 7r/8, 

£( c + 1 )]C]4p + 1 ~ c+1 - ( 5 - 5 °) 



1=0 



We substitute (jS33), do^oll and (l5~5TH) into (l5~33l to get that \f'/f\ < 
c + 1 + \n\ on dO z . Since c is equal to \k z \ or + 1 and \n\ < 3 (cf 
(|5.30jl ). the required estimate (|5.42|) follows. 



Step 4- Estimate of E^ x (u) and E^ y (u) We establish (|5.2fij) for j = x 
and j = y. For definiteness, we consider E^ x (u), and show that 

%M<1 + N 2 (5.51) 

(the calculations for En\ y {y) are essentially the same). On dO x , F = f, f is 
real and, from (j5.17|) . dw x /da = \ sec 2 a/2 < 1 for < a < n/2. It will be 
convenient to parameterise dO x by < w < 1 rather than by a. (j5.17j) and 
(f5~27Ijl give that 



f L f' 2 dw f 1 f' z 

= * / ^ ^ *. < 4 1 ^ da . (5.52) 

The estimate (j5.51j) requires more calculation than the corresponding 
result for j = z. It turns out that a pointwise bound on the integrand in 
()5.52|) is not sufficient, as f' 2 /(l + f 2 ) 2 3> \k x \ 3 on dO x . The domain on 
which f' 2 /{l + / 2 ) 2 >> |A; 2 | 3 has measure of order l/lfc^l, in keeping with 
(j5.51|) . But it will be necessary to estimate the integral in ()5.52|) itself. We 
note in passing that the complex representation (|5.15j) does not incorporate 
the cubic symmetries of the octant O in a simple way. Projecting along the 
axis (1,1,1)/ -\/3 rather than z would treat the boundaries symmetrically, and 
might lead to a simplification of the calculations below. 

To proceed, we collect the zeros and poles of / along dO x into a factor q 
(see (|5.filj) below for its explicit expression), writing 

/ = PQ, (5.53) 

where 

P = ABC = xfl(^±^y BC (5.54) 



/2 
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has no zeros or poles for < w < 1. We have that 







1 + p 





p 



+ p- 



p(pq + l/(pq)) l+p 2 q 2 
since \pq\ + l/\pq\ > 2, and 



< 



+max ( \p\ 



\P\ 



q 



1 + p 2 q 2 > 



l + q 2 , \p\>l, 
P 2 (l + q 2 ), \p\<L 



l + q 2 
(5.55) 



(5.56) 



With calculations similar to those in Step 3 (details are omitted), one shows 
that | (logp)'l is bounded on dO x independently of kj, ie 



< 1, < w < 1. 



(5.57) 



From (j5~22jl and JEH, if a = then p(0) = ±1, ie log |p(0)| = 0. If a ^ 0, 
using ()5.37|) we get that 



a/2 



|(logb(0)|)| = £ 



j=i 



log 



T2J-1 



a/2 



T2J 

a/2 



£>g(l-2/(a + 47))| 



j=i 



u / z /-2a 

< V4/(a + 4J)< / dx/(a + x) < In 3, (5.58) 
j=i Jo 



so that, in general, | (log |p(0)|)| < 1. Then (|5.57|) and ()5.58|) imply that, for 
< w < 1, | (log |pW|)| < I log(b(0)|)| + \p'/p\ dw' < 1, or 



max(|p|,l/H) <1, 0<w<l. 
Substituting ()5.59|) and ()5.57|) into ()5.55|) . we get that 



< 1 



l + q 2 



(5.59) 



(5.60) 



The explicit form of q is obtained from (j5.21j) . (j5.22aj) and ()5.53|) . 



q = w 



w — r 



(5-61) 



We partition 90^ into intervals whose endpoints are the zeros and poles of 
q and estimate q'/(l + g 2 ) on each. We consider in detail the interval Ij = 
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(r 2 j_i,r 2 j); the other intervals are treated similarly. Let x = 2a(w — r 2 j-\, 
so that x varies between and 1 on Ij. On Ij, we write 



q(w(x)) = w n (x) (gj(x)h(x)) 



P2J-1 



where 



h{x) 



x 



(5.62) 



(5.63) 



contains the zero and pole at the endpoints of Ij (an explicit expression for 
gj(x) is given in (|5.7()|) below). It is straightforward to show (the calculation 
is similar to that in ()5.55|) ) that 



1+q 2 



(w{x)) 



< l\n\ + 



( dw 
\ dx 



dgj/dx 



gj 



+ 



gjdh/dx 



1 + gjh 2 



(5.64) 



Substituting ()5.64j) into ()5.6()j) and noting that \n\ < 3, dw/dx = l/(2a) and 
dh/dx = —1/(1 — x) 2 , we get that 



<l + a 2 



dgj/dx 



fj.i 



(5.65) 



((i-^VI^I + x^jI) 2 ' 

The integral of the last term in (|5.65|) may be estimated as follows: 

/ ((1 - x) 2 /\gj\ + x 2 \gj\y 2 dx<Gj + Q~ 3 \ (5.66) 
Jo 



where 



mm gj , 

xe[o,i] 



Gj = max \gj\. 

xe[o,i] 



(5.67) 



To get (|5.66j) . suppose first that Qj < 1 < Gj. We divide [0, 1] into the three 
subintervals 



K=%\G-j\ L — [\G~j~ 1 , 1 — \Qj], M= 



(5.68) 



On K, the integrand in ()5.66|) is bounded by G 2 (l — x)~ 4 , so that the con- 
tribution from K to the integral is dominated by Gj. On L, the integrand 
is bounded by (2(1 — x)x)~ 2 (since, in general, \a/g\ + \bg\ > 2|a6| 1 / 2 ), so 
the contribution from L is dominated by Gj + Gj 1 - On M, the integrand is 
bounded by g^ 2 x~ A , so its contribution is dominated by Qj 1 . In case Gj < 1, 
let K = [0, |] and L = [f , 1 - Qj/2); in case Qj > 1, let L = [^Gj\ |] and 
M = [|,l]. 

From ()5.65|) and ()5.66|) it follows that 



1-2J 



/2 



T-2J 



/ 

-a (! + / 



dw < 1 + a max 



dgj/dx 



9j 



+ a(Gj + gj 1 ). 



(5.69) 
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We need to estimate the terms involving gj. From (j5.61|) and (|5.62j) . gj is 
given by 

gj(x) = Pj^(x)P a/ ^j(l - x), (5.70) 

where 



P N (x) 



N 

n 

K=\ 



x + 2K-i t(n + \x + \) ra x + 1; 



x + 2K T(N + I x +l)r(|x + |)' 



(5.71; 



Then leading-order asymptotics for T(z), ie \ogT(z) ~ (z — |) log 2; — 2; and 
(logr)'(z) ~ In 2; (see, eg, [I]), yields 



^iv(z) ^ and 

which in turn imply the estimates 
dgj/dx 



P' 



(x)<l for < x < 1, 



(5.72) 



< 1, < x < 1, 



9j 



((a/2 + 1 - J)/J) 1/2 , 1 < J < a/4, 

(J/ (a/2 + 1 - J)) 1/2 , a/4 < J < a/2. 

Substituting ()5.73|) into ()5.69|) . we get that 

f' 2 , , Ja/2 + 1 - 1/2 



2A2 



dw < 1 + a + a 



J 



+ a 



J 



a/2 + 1 - J 



(5.73) 



1/2 



(5.74) 



Estimating the contribution from the interval [r2j,r 2 j+i] to the integral 
(J5.65)) is carried out in much the same way. The differences are that i) h(x) 
is replaced by l/h(x) = —h(l — x) (which may be accommodated by the 
substitution x —> 1 — x) and ii) gj(x) in ()5.70|) is replaced by 



x + 2J - 1 



Pj-lWPafr-J-lO- 



X 



(5.75) 



x + 2 J — a 

But the expression in (|5.75|) and its logarithmic derivative satisfy the same 
bounds as do gj and dgj/dx in ()5.73j) . Thus, the integral of / /2 /(l + / 2 ) 2 
over [r 2 j,r 2 j+i] satisfies the same bound ()5.74j) as does the integral over 
[ r 2j-i ; r 2j]- We obtain a bound on the collective contribution from the inter- 
vals [rj, by summing over J in ()5.74|) . 



/2 



a/2 



2^2 



dw < l+a 2 +a 2^ 



j=i 



a/2 + 1- J 
J 



1/2 a/2 



1/2 



^ Va/2 + 1 - J / 

(5.76) 
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The first sum may be estimated as 

(5.77) 

and the second is similarly bounded. Thus 



f 



r a f/2 



f 

dw<l + a 2 . (5.78) 



;i + / 2 ) 2 



The contributions from the remaining intervals [0, r%] and [r a , 1] are treated 
similarly, and we get 

f' 2 

dw<l + a 2 , (5.79) 



;i + /2)2 
f 

2\2 



1 + / 



dw < 1 + a 2 . (5.80) 



We give an argument for ()5.79|) ( ([5.80)1 is treated similarly). For definiteness, 
let us assume that pi = 1 (the case pi = — 1 is treated similarly). On [0, r±\ 
we write q = uv, where v = w n (w — r\) contains the zeros at the endpoints 
and 

it = {w r^p-^l - 2a(w - f )) (5.81) 

contains the remaining factors. Arguing as in (J5.64)) and 1)5.60)1 (but noting 
that u and v are functions of w, not a rescaled coordinate x), we get that 

f 2 ^ u' 2 u 2 v' 2 , . 

£ 1 + -T + 77-—^- (5-82) 



(l + /2)2 ~ M 2 ( 1+M 2 W 2)2' 

Clearly v /2 < 1, so that u 2 v' 2 /(l + u 2 v 2 ) 2 < u 2 . From ()5~72)l and JEHU), 
|u'/u| 2 < a 2 and |u| 2 < a. Therefore, f 2 /(l + / 2 ) 2 < 1 + a 2 , and l|On|l 
follows. The required bound on E^ x (u), 1)5.51)1 . follows from substituting 
the estimates (I5~T%1) - (l5~£U)l into the formula (15321) . 

□ 

Remark. To extend the upper bound of Theorem El to, say, a general convex 
polyhedron, one would like to have a generalisation of the octant configura- 
tions of Section 15.21 These would be conformal maps, perhaps with singu- 
larities, of a general convex geodesic polygon £ C S 2 into S 2 such that each 
edge of £ is mapped into the geodesic which contains it. 
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A Wrapping numbers as complete invariants 

In [21] we gave a homotopy classification of tangent unit-vector fields on 
a convex polyhedron P in terms of a set of invariants called edges signs, 
kink numbers and trapped areas. Here we show that these invariants can be 
determined from the wrapping numbers, so that the wrapping numbers w aa 
constitute a complete set of invariants. 

It suffices to consider the invariants associated with a single vertex. Let 
/ be the number of faces of P. Let v denote a vertex of P, and let w a denote 
the wrapping numbers on a cleaved surface around v (we suppress the vertex 
label a), where a is an /-tuple of signs. Suppose v has b > 3 coincident 
faces and therefore b coincident edges. Let E r , 1 < r < b, denote the edges 
coincident at v, ordered consecutively clockwise with respect to a ray from 
v through the interior of P. By convention let E b+1 = E 1 . Let E r denote 
the unit vector along E r directed away from v. Let F r denote the face with 
edges E r and E r+1 . An (unnormalised) outward normal on F r is given by 

F r = x gr^ (A.l) 

Let us explain briefly how the edge signs, kink numbers and trapped areas 
are defined (see [21] for details). Fix a homotopy class C(h), and let n e C(h) 
denote a representative. The edge sign e r is given by the orientation of n 
along the edge E r relative to E r , so that n(r) = e r W for r G E r . The kink 
number k r is an integer giving the winding number of n about F r along a 
path on the face F r starting on the edge E r+1 and ending on the edge E r 
(tangent boundary conditions imply that n is orthogonal to F r along such 
path). A minimal winding (eg, n taking values along the shortest arc from 
E r+1 to E r on S 2 ) has kink number equal to zero. Finally, the trapped area 
Q is the area in S 2 of the image under n of a cleaved surface about v (the 
area is normalised so that S 2 has area An). 

It is straightforward to derive the following expression for the wrapping 
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numbers in terms of e r , k r and Q |21j : 

b 

Q = 47rw a(s) + 2tt s g n ( Fr • S ) F 

r=l 

6-1 

+ ^ (A^E 1 , e r E r , e r+1 E r+1 ) - 4ttt(s; e^ 1 , e r E r , e r+1 E r+1 )) . (A.2) 

r=2 

Here, s 6 5 2 may be taken to be any unit vector which is transverse to every 
face of P (including faces which are not coincident at v). er(s) gives the sector 
to which s belongs. For a, b, c G S 2 , the quantities A (a, b, c) and r(s; a, b, c) 
are defined as follows. Let K G S 2 denote the spherical triangle with vertices 
a , b, and c (K is well defined provided a, b and c are not coplanar and no 
pair of them are antipodal). Then A(a, b, c) is the oriented area of K with 
values between — tt and tc (the sign is given by sgn (a-(bxc))). The quantity 
r(s; a, b, c) is given by 

i u \ f s g n ( a - ( b x c )). seK, 
>a '"' j„. , iK . (A ' 3) 

That is, r(s; a, b, c) is equal to zero unless s belongs to K, in which case it 
is equal to ±1 according to whether K has positive or negative area. Note 
that r(s; a, b, c) is well defined if s is transverse to the planes spanned by 
a, b, c taken pairwise. 

Our task here is to show that, given the wrapping numbers w a for topol- 
ogy h, we can determine the edge signs, kink numbers and trapped areas. 
We begin by determining the edge signs, specifically e r and e r+1 . Without 
loss of generality, we can assume that r / 1 and r ^ b (note that (jA.2|) 
remains valid if the edge indices are cyclically permuted; there is nothing 
special about the edge E 1 ). 

Let S C S 2 denote the great circle containing E r and E r+1 . The four 
points ±E r , ±E r+1 partition S into four disjoint open arcs. Denote these 
by S m , m = 1,2,3,4 (the ordering is not important). Let S mt denote the 
arc whose endpoints are e r E r and e r+1 E r+1 . We determine m*, and hence 
e r and e r+1 , by means of the following calculation. For each m, choose some 
s m € S m , and let 

s± = s m ± eF r (A.4) 

(recall that F r is normal to S) with e > small enough so that s^j is trans- 
verse to every face of P and so that 

sgn (F s ■ s+) = sgn (F s ■ s m ), s + r. (A.5) 
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We subtract the two equations obtained by letting s = s^j in (jA.2|) to 
obtain 

- k r + r(s+ ; e^ 1 , e r E r , e r+1 E r+1 ) - r(s m ; e^ 1 , e r E r , e r+1 E r+1 ). (A.6) 

Let K C S 2 denote the spherical triangle with vertices e 1 E 1 , e r E r and 
e' r+1 E r+1 . If to ^ m*, then neither s+ nor s~ lies in K, so that, by (|A.3|) . 
both of the r-terms in (|A.6|) vanish. On the other hand, if to = m*, then 
either sj£ or s~ lies in K but not both, so that one of the r-terms in flA.6|) 
vanishes while the other is equal to ±1. Therefore, amongst the four pos- 
sible values of w ff ' Sm ' — w a ^ Sm \ three will have the same value and one will 
be different, m* is identified as the index for which w' J ' Sm ' — w a ^ Sm ^ has the 
different value. 

Once the edge signs are determined, the kink numbers can be obtained 
from (|A.6j) . and hence the trapped from (IA.2I). 
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